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We apply the conformal gravity theory to a sample of 111 spiral galaxies whose rotation curve 
data points extend well beyond the optical disk. With no free parameters other than galactic mass 



> 

in 

■ to light ratios, the theory is able to account for the systematics that is observed in this entire set of 

cn: 

rotation curves without the need for any dark matter at all. In previous applications of the theory 
a central role was played by a universal linear potential term V{r) = that is generated 

through the effect of cosmology on individual galaxies, with the coefficient 70 = 3.06 x 10^'^'^ cm^^ 
^ , being of cosmological magnitude. Because the current sample is so big and encompasses some 

I specihc galaxies whose data points go out to quite substantial distances from galactic centers, 

we are able to identify an additional globally induced universal term in the data, a quadratic 
V{r) = —K(?r'^ jl term that is induced by inhomogeneities in the cosmic background. With k 
being found to be of magnitude k = 9.54 x 10~^^ cm~^, through study of the motions of particles 
contained within galaxies we are thus able to both detect the presence of a global de Sitter-like 
component and provide a specihc value for its strength. Our study suggests that invoking dark 
matter may be nothing more than an attempt to describe global physics effects such as these in 
purely local galactic terms. 
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I. INTRODUCTION 



Observational studies of spiral galaxies have repeatedly established that galactic rota- 
tional velocities look nothing like the velocities that would be produced by the Newtonian 
gravitational potentials associated with the luminous matter in the galaxies. In consequence, 
it is quite widely thought that such velocity discrepancies are to be explained by the presence 
of copious amounts of non-luminous or dark matter in galaxies. Since the case for the pres- 
ence of such dark matter rests solely on the assumption that wisdom acquired from studies 
on solar system distance scales can be extrapolated without modification to the much larger 
galactic distance scales, a few authors have ventured to suggest (see e.g. [1] for a recent 
review) that dark matter may not actually exist and that instead it is the standard New- 
tonian description that needs modifying. In this work we apply one particular candidate 
alternative theory, namely conformal gravity, to a large and comprehensive sample of 111 
galactic rotation curves. With only one free parameter per galaxy, the galactic mass to light 
ratio, we find that the conformal theory provides for a good accounting of the data without 
the need for any dark matter at all. Moreover, because our sample is so large, through our 
fitting we are able to find evidence in the data for the presence of a universal quadratic 
potential term that the conformal theory possesses. 

As a theory, conformal gravity (see e.g. |l|) is a completely covariant metric theory 
of gravity that possesses all the general coordinate invariance and equivalence principle 
structure of standard Einstein gravity, but which in addition possesses a local conformal 
invariance in which the action is left invariant under local metric transformations of the 
form gfj,u{x) — ?■ e^"'^^^ g ^^{x) with any arbitrary local phase a{x). As a symmetry, confor- 
mal invariance forbids the presence of any fundamental cosmological constant term in the 
gravitational action, with the action being uniquely prescribed by the Weyl action 

where 

Cx^uK = RxflUK — 2 (gxuRfiK — gxuRfiu — gfiuRxn + gfinRxv) + g-^"a {gxugfiK — gxngiiu) (2) 

is the conformal Weyl tensor and the gravitational coupling constant a, is dimensionless. 
Thus, unlike the standard Einstein-Hilbert action /eh = — (I/IOttG) / d^x (—gY^'^R"^, which 



R,^R^- - {l/3)iR\y , (1) 
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can be augmented to include a / d'^x {—gY^'^A term, the conformal theory has a control over 
the cosmological constant that the standard Einstein theory does not, and because of this 
one is able to provide a potential solution to the cosmological constant problem 



II. LOCAL CONSIDERATIONS 



For the Weyl action the equations of motion take the form 



1^(2) 3 1^(1) 



T"", (3) 



where 



= 2g^^{R\y^.^p - 2{R\y^-^'' - 2R\R^^ + -g^^'iR^f, 

= \9"'{R\y^-B + R'^'^-B - R'^'^-B - R'^^'^-B - 2R^^R% + \g^'R^pR^^ (4) 



Thus, since W^" vanishes when R^" vanishes, we see that, as well as being a vacuum solu- 
tion to Einstein gravity, the Schwarzschild solution is also a vacuum solution to conformal 
gravity. The conformal theory thus recovers all the standard solar system Schwarzschild 
metric phenomenology, just as is needed for any metric theory of gravity. 

However since the vanishing of W^^ could potentially be achieved without R^'^ needing 
to vanish, the conformal theory could also have some non-Schwarzschild solutions as well. 
To determine what such solutions might look like, Mannheim and Kazanas solved for the 
metric outside of a localized static, spherically symmetric source of radius tq embedded in a 



region with T^i/(r > tq) = 0. They found j4| that in the conformal theory the exact, all-order 
classical line element is given by ds"^ = —B{r)c^dt^ + dr"^ / B{r) + r'^dQ2 where the exterior 
metric coefficient B{r > tq) is given by 

2/3 

B{r > ro) = w \- — kr^ (5) 

with w = (1 — 6/37)^''^. In equation (|5]) the presence of the three integration constants /3, 7 
and k is due to the fact that unlike the standard second-order derivative Einstein theory, the 
conformal theory is instead based on fourth-order derivative equations, to thus contain two 
additional terms. With the /37 product numerically being found to be small for standard 
astrophysical sources (see below), we can set w = 1. Then with the emergence of a 1 — 
2/3/r term we see that the conformal gravity metric contains the familiar general-relativistic 



Schwarzschild metric solution (and thus its non-relativistic Newtonian gravitational limit as 
well), while departing from it only at large r, i.e. departing from it in precisely the kinematic 
region where the dark matter problem is first encountered. 

In seeking to relate the various integration constants in equation to properties of 
the energy-momentum tensor T^^, of the source, Mannheim and Kazanas found [5] that in 
terms of the general source function /(r) = (3/4ac,i?(r))(T°Q — T\), the exact fourth-order 
equation of motion given in equation ([3]) reduced to the remarkably simple form 

-{W% - ly;) = V^B = B"" + ^ = = /(r), (6) 



B{ry ^ ^' • ■ r r 

without any approximation whatsoever. (The primes here denote derivatives with respect to 
r.) Since V^(r^) vanishes identically everywhere while V^(l/r) and V^(r) evaluate to delta 
functions and their derivatives, we see that of the integration constants given in equation 
([5]), only 13 and 7 can be associated with properties of a local source of radius tq; with the 
matching of the interior and exterior metrics yielding 

7 = r dr'r''f{r'), 2fi = \ T dr'r'^j{r'). (7) 

/Jo D JO 

Since the —kr^ term in equation ([5]) is a trivial vacuum solution (T^i,(r) = everywhere) 
to equation (E]), as such it is not coupled to the local matter source, with the discussion 
here providing no basis for considering it further. (As we will show below, once we allow 
for matter sources in the r > tq region (i.e. T^i/(r > tq) 7^ 0), we will be able to generate a 
quadratic term that will be coupled to the matter in the r > ro region.) 

With the /3 and 7 terms coupling to the local source, we see that in conformal gravity a 
given local gravitational source generates a gravitational potential 

V*{r) = -^ + l^ (8) 

per unit solar mass, with /3* being given by the familiar MqG/c^ = 1.48 x 10^ cm, and with 
the numerical value of the solar 7* needing to be determined by data fitting. For V*{r) to 
generate non-relativistic motions, it is necessary that /3*/r and 7*r both be very much less 
than one, and that the /3*7* product thus be very much less than one too. The domain in 
which we can use V*{r) is thus an intermediate one where r is neither too small nor too 
large, and in which /3* and 7* are such that /3*7* (and thus the quantity w ~ 1) are very 
much less than one. In our applications of V*{r) to galaxies we shall find that all of these 



conditions hold. However, since the form B{r > tq) = w — 2/3/r + 7r found above is exact 
without approximation, in the event that we were not in the domain needed to use V*{r), 
we would have to use the exact geodesies associated with the exact B{r) instead. 

In conformal gravity the visible local material in a given galaxy would generate a net 
local gravitational potential Vloc(^) given by integrating V*{r) over the visible galactic 
mass distribution. Typically, the luminous material in a disk galaxy is distributed with a 
surface brightness = Soe~^/^° with scale length Rq and total luminosity L = 2nT.oRl, 
with most of the surface brightness being concentrated in the R < ARq or so optical disk 
region. For a galactic mass to light ratio M/L, one can define the total number of solar 
mass units A^* in the galaxy via {M/L)L = M = N*Mq. Then, on integrating V*{r) over 
this visible matter distribution, one obtains jl| the net local luminous contribution 



vloc ^ N*f3*c'R 
R 2Rl 



/i?\ /-R\ / R \ / R 



^'■'■^\(^]kJ^) (9) 
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for the centripetal accelerations of particles in circular orbits in the plane of the galactic 
disk. In the R^ Rq limit this expression simplifies to 

vloc ^ N*P*c' (,,^\, N*Yc' A _ 3i?| _ A5Rt\ N*P*c' N*^^ 
R i?2 \^ 2R^)^ 2 \ 2i?2 ) R^ 2 ' ^ > 

with the entire galaxy acting as if it were a point source located at the galactic center. With 
the surface brightness of the optical disk region essentially becoming negligible hy R = 4i?o 
or so, equation ( ITOj) can be expected to be a good approximation to equation ([9]) at points 
with R > ARq. 

III. GLOBAL CONSIDERATIONS 

Unlike the situation that obtains in standard second-order gravity, one cannot simply use 
equation (Q as is to fit galactic rotation curve data, as one must take into consideration the 
effect of the rest of the material in the Universe as well. To see why this is the case, we recall 
that for standard gravity, the solution to the second-order Poisson equation V^0(r) = g{r) 
for a general static, spherically symmetric source g{r) is given by 

(j){r) = — j dr' r''^g{r') — J dr'r'g{r'), (11) 



with derivative 

"^^^'^ ^ t'dr'r-g{r'). (12) 



dr Jo 

As such, the import of equation (fT2|) is that even though g{r) could continue globally all the 
way to infinity, the force at any radial point r is determined only by the material in the local 
< r' < r region. In this sense Newtonian gravity is local, since to explain a gravitational 
effect in some local region one only needs to consider the material in that region. Thus in 
Newtonian gravity, if one wishes to explain the behavior of galactic rotation curves through 
the use of dark matter, one must locate the dark matter where the problem is and not 
elsewhere. Since the discrepancy problem in galaxies occurs primarily in the region beyond 
the optical disk, one must thus locate galactic dark matter in precisely the region in galaxies 
where there is little or no visible matter. 

Despite the fact that the force in equation (IT^ is not sensitive to any material beyond the 
radial point of interest, this local character to Newtonian gravity is not a generic property 
of any gravitational potential. In particular for the fourth-order Poisson equation V^0(r) = 
h{r) = /(r)c^/2 of interest to conformal gravity, the general solution is of the form 

0(r) = — / dr'r'^h{r') / dr' r'^hir') - - / dr' r'^hir') / dr' r'hir'). (13) 
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With the derivative of the potential evaluating to 

= -1 r dr'r'^h{r') + ^ [ dr' r'%{r') - - dr' r'h{r'), (14) 
dr 2, Jo 6r j o 3 Jr 

this time we do find a contribution to the force coming from material that is beyond the 
radial point of interest. Thus in the third integral in equation ( fT4l) we recognize a potential 
global contribution to local motions, with a test particle in orbit in a galaxy being able to 
sample both the local field due to the matter in the galaxy and the global field due to the 
material in the rest of the Universe as well. In conformal gravity then, to determine motions 
of particles inside of galaxies one cannot ignore the effect of the material outside of them. 

In order to determine the effect that material exterior to galaxies might have on galaxies, 
we note that there are actually two global effects that we need to take into consideration. 
Specifically, we need to consider the effects of both the homogeneous background cosmology 
and the inhomogeneities that are present in it. Moreover, in the conformal theory, these ef- 
fects have very different geometric structures. The global background cosmology is described 
by a comoving Robertson- Walker (RW) geometry. Since an RW geometry is homogeneous 



and isotropic, its metric is conformal to flat. Thus in an RW geometry both the Weyl ten- 
sor and W^'^ vanish identically. However, since by their very nature inhomogeneities can 
localize in space, they are associated with geometries in which neither the Weyl tensor nor 
W^'^ can vanish. Indeed, in the derivation of equation given in j5|] , it was found that in a 
static, spherically symmetric geometry the quantity (3/i?(r)) {W^q — W^^) evaluates exactly 
to \/^B{r), with some components of W^'^ necessarily being non-zero in any conflguration 
in which V^-B(r) is non-zero. Thus it is only inhomogeneities that contribute to the third 
integral in equation (fT^ . 

As regards the cosmological background, we note that since the background is associated 
with W^" = 0, it too will contribute to the solution to V^-B(r) = /(r). However, it will 
do so not as part of the particular integral solution given in equation f[T^ but as part of 
the complementary solution to V^-B = instead. The background cosmology can thus have 
physical consequences for galactic motions provided it causes W'^'^ to vanish non-trivially, 
i.e. provided it causes T^'^ to vanish non-trivially in equation ([3]). We thus recall a] 
that in conformal gravity one can indeed construct cosmologies in which T'^'^ does vanish 
non-trivially, and in them the scale factor R{t) and the 3-curvature K of the RW metric are 
related to the cosmological matter content, with K being found [Ij to be negative. 

With the Hubble flow being described in comoving coordinates and galactic rotational 
velocities being measured in a coordinate system in which a galaxy is at rest, to determine the 
effect of the Hubble flow on galactic motions we need to transform the RW metric to a static 
coordinate system. To this end, we recall j] that the general coordinate transformation 

4r 



P 



2(l + 7or- A;r2)V2 + 2 + 7or 
effects the metric transformation 



T 



dt R{t) 



(15) 



-(1 + 7or - kr^)c^dr + 



(1 + 7or — kr"^) 



+ r^dn2 



1 [1 - 7oV/16 - A;pV4]2 
R^{t) [(l-7op/4)2 + A;pV4]^ 



-c^dr' 



R\r) 



[l-(7oVl6 + fc/4)p^ 



212 



(dp" + p'dVt^ 



With the transformed metric being written compactly as 



ds' 



,2a(T,p) 



-c'dr' + 



RH 



dp^ + p^dQ.2 



.(16) 



(17) 



[1 + Kp'^/AY 

we see that the transformed metric is conformally equivalent to a comoving RW metric as 
written in spatially isotropic coordinates with spatial 3-curvature K = — 7o/4 — k. Since an 



RW geometry is conformal to flat and since it remains so under a conformal transformation, 
we see that when written in a static coordinate system a comoving conformal cosmology 
looks just like a static metric with universal linear and quadratic terms. 

With the dynamics that leads to the K < RW metric in the first place only involving 
one physical cosmological scale and not two [1], the decomposition of just one RW scale (viz. 
K) into two static scales (70 and k) is artificial, and it is not meaningful to keep both 70 
and k. With a K < RW metric only being defined for p < 2/(— fC)^/^, we see that if we 
keep the k term alone, which would require k = —K to be positive, the coordinate r would 
become complex at r = with a transformation of the form p = 2r/[(l — kr'^Y^'^ + 1], 

r = p/(l + kp^/A) only being able to cover an r < k^^^"^ patch of a static spacetime 
geometry. However, if we keep the 70 term alone in equation f[T^ . then provided 70 is taken 
to be positive, the coordinate r would then be able to run all the way to infinity (just as one 
would want of a static, spherically symmetric geometry), with nothing being able to become 
complex. Thus, to be able to recover the standard static geometry with its infinite domain 
for the coordinate r, we shall retain the 70 term and leave out the k term in the coordinate 
transformation. (For a discussion of a possible role for the /c-dependent transformation in 
a cosmological context see [7].) And while we shall uncover yet another quadratic term 
below, viz. the — /tr^ term that is associated with cosmological inhomogeneities, we will find 
that its scale is sub-cosmological and thus not to be associated with the cosmological — fcr^ 
term given in equation (fT6!) . (For purely phenomenological purposes, if one were to keep 
both of the cosmological 7or and —kr"^ terms in equation ( fT6i) anyway, one could consider the 
quadratic term to be used in the fitting described below to be a composite of homogeneously 
and inhomogeneously induced quadratic terms.) 

On dropping k we replace equations f|T5l) and f|T6l) by 

4r p 



P 



and 



2(l + 7or)i/2 + 2 + 7or' (1 - 7op/4)2 



;i + 'yor^dt^ + — + r^dQ^ 

(1 + lor) 



J dtR{t) (18) 



1 /l + 7op/4\ 



i?2(r) \ I-70P/4/ 



(19) 



[l-7oV/16]^ 

Without the k term the RW 3-curvature is given hj K = — 7q/4, a necessarily negative quan- 
tity. Since the only way to make K be positive would be to have complex 70, and the only 
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way to make K be zero would be to have 70 = 0, we see that in the rest frame of a comoving 
galaxy (i.e. one with no peculiar velocity with respect to the Hubble flow), a topologically 
open comoving cosmology (viz. just the one found in and only a topologically open 
one, looks just like a universal linear potential, with a strength given by 7o/2 = (—KY^'^. 

In the conformal theory then we recognize not one but two linear potential terms, a local 
]V*7*-dependent one associated with the matter within a galaxy and a global cosmological 
one 7oc^r/2 associated with the cosmological background. Thus in [s^ it was noted that in 
the weak gravity limit one could add the two potentials, with the total circular velocity f tot 
then being given by 

Vtot = ^Loc + ' (20) 

with asymptotic limit 

Vtot ^ + ^ + • (21) 

In equation (!20|) was used to fit the galactic rotation curve data of a sample of 11 galaxies, 
and good fits were found, with the two universal linear potential parameters being found to 
be given by 

7* = 5.42 X 10"^^ cm~\ 70 = 3.06 x 10"^° cm'^ (22) 

The value obtained for 7* entails that the linear potential of the Sun is so small that there 
are no modifications to standard solar system phenomenology, with the values obtained for 
A^*7* and 70 being so small that one has to go all the way to galactic systems before their 
effects can become as big as the Newtonian contribution. Moreover, the value obtained for 
7o shows that it is indeed of cosmological magnitude, just as desired. 

While the analysis described above provides no unequivocal reason for including any 
possible quadratic potential term in equation fl20|) . valid justification for considering it is 
obtained by considering not the homogeneous cosmological background, but rather the in- 
homogeneities in it. On large scales these inhomogeneities would typically be in the form 
of clusters and superclusters and would be associated with distance scales between 1 Mpc 
and 100 Mpc or so. Without knowing anything other than that about them, we see from 
equation ( JT3l) that for calculating potentials at galactic distance scales (viz. scales much 
less than cluster scales) the inhomogeneities would contribute constant and quadratic terms 
multiplied by integrals that are evaluated between fixed end points, to thus be constants. 
(I.e. all that we require of the — (r^/6) dr' r'h{r') integral in equation (fT3!) is that it begin 
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at some minimum cluster-sized radius rdus that is outside the galaxy and independent of it.) 
Thus given the quadratic term in equation f|T3|) . then again up to peculiar velocity effects, 
for weak gravity, and on scales r < rdus, "we can augment equation ( 12U]) to 

%OT = ^Loc + — i? (23) 

where k(? = (1/3) J^^^^dr' r'h{r'), with associated asymptotic limit 

As such, equation ( 123|) can be derived from a metric with a term B{r) = — /tr^, and thus has 
a de Sitter-like form. However, it is not associated with an explicit de Sitter geometry per se 
since the inhomogeneities that give rise to it are not distributed in a maximally 4-symmetric 
way. Nonetheless, a particle in orbit in a galaxy would be affected by the quadratic term, and 
thus behave in exactly the same way as if it had been embedded in a de Sitter background. 

Now in a conformal theory particles can only acquire mass through some scalar field 
symmetry breaking procedure, and thus when particles propagate they can exchange energy 
and momentum with such fields. However, as we show in the Appendix, even in the presence 
of such an exchange, circular orbits in galaxies are still of the geodesic form rd(j)/dt = 
[rc^-B'(r)/2]^/^ that leads to equation ( |23|) . Equation (!23|) with its universal k is thus our 
main theoretical result, and so we proceed now to apply it to galactic rotation curve data. 



IV. CONFORMAL GRAVITY DATA FITTING 

Since successful rotation curve fitting to an 11 galaxy sample was obtained in 8] via 
the use of equation (l2Up . one would initially anticipate that even if the —kc^R term in 
equation (l23l) were to be present in principle, in practice it would be too small to have any 
effect. However, the sample we study here is much larger (111 galaxies) and it contains 
some galaxies whose data points extend to far larger distances 
had been the case for the 11 galaxy sample originally studied in 
through fitting 21 such highly extended galaxies that we were able to uncover a role for the 
—KC^R term and extract a value for k given by /t = 9.54 x 10~^^ cm~^ ^ (100 Mpc)~^. In the 
fitting to the full 111 galaxy sample we shall use this value for k and the values for 7* and 
7o as given above in equation (122|) . For the fitting then there is just one free parameter per 
galaxy, namely the galactic mass to light ratio, and thus our fitting is highly constrained. 

10 



Tom galactic centers than 
8|. As reported in [9|, it is 



With the stars in galaxies lying within the optical disk region, to fully explore the rotation 
curves of galaxies one needs to study the HI gas spectra as it is only the gas in galaxies that 
extends well beyond the optical disk region. To get velocity measurements that are free of 
projection concerns one wants galaxies to be close to edge on along our line of sight, and 
to be able to model the gravitational contribution of the luminous disk one needs good disk 
photometry. Given these criteria there is a now quite substantial number of galaxies for 
which one can do modeling, with the 111 galaxy set that we use being a large, very varied 
and representative sample that contains both high surface brightness (HSB) galaxies where 
both A^* and Sq are large, and low surface brightness galaxies with small Sq and dwarf 
galaxies with small A^* (collectively referred to here as LSB galaxies since many small N* 
galaxies have small Sq and vice versa). 

Having this broad a variety of galaxies turns out to be very instructive since one of the 
most interesting aspects of equations ( 120!) and (123|) is that there are situations in which 
departures from the luminous Newtonian prediction can be very pronounced. One situation 
is when A^* is small, since then the net Newtonian contribution cannot compete with the 
fixed magnitude 70 and k terms. Another situation is when the quantity N*/Rq ~ Sq is 
small. Specifically, since the Newtonian contribution in equation (Q (the /3* dependent 
IqKq — IiKi term) numerically peaks at around R = 2.2Ro, the strength of the Newtonian 
term at the peak will be set by the magnitude of N*/Rq, and when small will not be able 
to compete with the fixed magnitude 70 and k terms. Since the linear term dominates over 
the quadratic one until the largest distances, in both small A^* and/or small So galaxies one 
should expect the rotation curves to start rising immediately, just as is systematically seen 
in the data sample. The case where the luminous Newtonian contribution is not suppressed 
is in HSB galaxies, and here the falling Newtonian contribution can compete with the rising 
linear term to give a region of approximate flatness before any rise could set in, again just 
as is systematically seen in the data. Thus we see that the simple formula given in equation 
(1231) directly captures the essence of the data, and as the fits show, the formula captures not 
just the qualitative trend but the actual quantitative numerical values of the velocities as 
well. Finally, we note that for all galaxies the quadratic term will eventually take over, to 
then arrest the rising linear potential terms and cause all rotation velocities to ultimately 
fall. Moreover, since v"^ cannot go negative, beyond R ~ {N*'~f* + 7o)/2k (~ 100 kpc for 
A^* = 70/7* = 5.65 * 10^°) there could no longer be any bound circular orbits, with galaxies 
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thus having a natural way of terminating, and with the allowable sizes of galaxies being 
determined by an interplay between galaxies and the global structure of the Universe. 

For the actual fitting we have predominantly used galaxies that were studied in large 
surveys. In particular for the rotation curves we have used 18 galaxies from THINGS: The 
HI Nearby Galaxy Survey (as detailed in Table 1), 30 galaxies from a study of the Ursa Major 
Cluster of Galaxies (Table 2), 20 galaxies from a study of LSB galaxies, as augmented by 
an extended distance study of UGC 128 (Table 3), 21 galaxies from a second study of LSB 
galaxies (Table 4), and also included some 22 miscellaneous galaxies (Table 5), with this last 
set containing many of the galaxies that played a significant historical role in establishing 
that there actually was a galactic missing mass problem in the first place. The sample we 
use contains all the 11 galaxies that were studied in (DDO 154, DDO 170, NGC 1560, 
NGC 3109, UGC 2259, NGC 6503, NGC 2403, NGC 3198, NGC 2903, NGC 7331, and 
NGC 2841), with a few of them having undergone significant updates since then. Of the 
lllgalaxies in our sample, the 21 that extend the furthest in radial distance were reported 
in 9|, and for completeness we also include them here. In order of increasing largest radial 
distance the 21 galaxies are NGC 3726, NGC 3769, NGC 4013, NGC 3521, NGC 2683, UGC 
1230, NGC 3198, NGC 5371, NGC 2998, NGC 5055, NGC 5033, NGC 801, NGC 5907, NGC 
3992, NGC 2841, UGC 128, NGC 5533, NGC 6674, UGC 6614, UGC 2885 and Malin 1. 

For the fits we have taken photometric luminosities, optical disk scale lengths and HI gas 
masses from Refs. 10| through 66j. The values we use are listed in Tables 1-5. In the last 
column in each of these Tables each set of four references gives the data sources for rotation 
velocities (f ), luminosities (L), disk scale lengths (i?o) and HI gas mass (HI). 

As described in the Appendix, for 11 of the galaxies (NGC 801, NGC 2998, NGC 5033, 
NGC 5055, NGC 5371, NGC 5533, NGC 5907, NGC 6674, UGC 2885, ESQ 1440040 and 
Malin 1), we have also included the contribution of a central spherical bulge. Since HI gas 
distributions extend well beyond optical disk distributions, for simplicity we have modeled 
the gas profile in each galaxy as a single exponential disk with a scale length larger than that 
of its optical disk. And on finding little sensitivity to the actual gas scale lengths used (since 
the gas makes only a small contribution to rotation velocities), for definiteness we took the 
ratio of gas scale length to optical disk scale length in each galaxy to be four. Also we have 
multiplied the overall HI gas contribution by 1.4 to account for primordial Helium. (When 
an HI gas mass was not available, the HI gas mass is listed as NA in the Tables.) In the fits 
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the gas contribution is never that significant. Specifically, in the HSB galaxies the mass in 
stars is much greater than the mass in gas, while in the LSB galaxies, neither the gas nor 



the stars are ab 
discussion in 



e to compete with the universal 70 and k terms. In the fits we followed the 
55| and required that M/L not be less than 0.2Mq/Lq. In the Tables we have 
listed the fitted stellar mass to light ratios {M/L) that we have obtained from our fitting, 
with the (M/L)stars values quoted in the Tables representing the total stellar disk plus bulge 
mass combined as divided by the total blue galactic luminosity in those cases where we have 
included a galactic bulge. In almost all cases the mass to light ratios that we obtain are 
reasonably close to the mass to light ratio found in the local solar neighborhood, just as one 
would want. 

In those cases where optical scale lengths have been measured in many wavelengths, by 
and large we have used the scale lengths as measured in the longest available wavelength 
band (usually the K band) and have systematically done so for the entire 30 galaxy Ursa 
Major sample. For seven of the galaxies (NGC 7137, UGC 477, ESQ 840411, ESQ 1200211, 
ESQ 3020120, ESQ 3050090 and ESQ 4880490) little or no surface photometry is available 
at all. As described in the Appendix, for these particular galaxies we have had to estimate 
scale lengths, with the sources for the scale lengths for these galaxies accordingly being listed 
as ES in the Tables. 

The place in our theory where there is the most sensitivity to parameters is in the adopted 
distances to the individual galaxies, since the parameters 7*, 70 and k that appear in equation 
( 123|) are given as absolute quantities. To establish a common baseline for determining 
adopted distances, for all the galaxies in our sample we have used the distances listed in the 
NASA/IPAC Extragalactic Database (NED). In this database distances are obtained either 
via direct visual measurements (typically Cepheids or the Tully-Fisher relation) or indirectly 
via redshift measurements. For the directly determined distances a world average mean value 
and its one standard deviation uncertainty are listed. The redshift-based determinations 
depend on how one models both the peculiar velocity with respect to the Hubble fiow of the 
Milky Way Galaxy and the peculiar velocity of the galaxy of interest. Five different such 
models are provided in the NED, and with each one giving a mean value and uncertainty, 
taken together the five determinations and their uncertainties provide a spread in values. 
For definitiveness, for redshift-based distance determinations we have opted to use the mean 
value associated with the galactocentric distance determination. For our entire set of 111 
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galaxies there was only a handful of 10 galaxies for which using the visually-determined mean 
or the redshift-determined galactocentric mean did not immediately give a reasonable fit. 
For IC 2574, NGC 2403, NGC 3621, NGC 7793 and NGC 3109 we found it advantageous 
to use adopted distances up to one standard deviation above the NED mean, while for 
NGC 2841, DDO 170, NGC 5033 and NGC 5533 we allowed up to one standard deviation 
below the NED mean. (For NGC 2841 the adopted distance we used coincides with the one 
given by Cepheid data alone, with TuUy-Fisher based determinations yielding a somewhat 
higher value.) For NGC 6674 we used the smallest allowed distance value within the redshift 
determined spread in values. Thus for no less than 101 of the galaxies in our sample our 
theory captured the essence of the rotation curve data using the NED preferred distances 
as is. And moreover, despite the fact that the sensitivity to adopted distance is the most 
pronounced in the 21 large galaxy sample, for only four of them (NGC 5033, NGC 2841, 
NGC 5533 and NGC 6674) did we even need to consider not using the NED mean values as 
is. The fact that our fits work so well at the NED distances is thus a noteworthy achievement 
for our theory. In the Tables we have listed the specific adopted distances that we have used. 

Of the galaxies we fit, the data for Malin 1, a giant LSB galaxy, go out further in radial 
distance than any other of the galaxies in our sample, and as such this galaxy actually 
provides the sternest test of our ideas. Malin 1 is unusual in that its adopted distance is 
far larger than that of any other galaxy in our sample. In fact it is at such a high redshift 
{z = 0.0824) that its luminosity and angular diameter distances differ quite significantly. 
For Malin 1 the NED gives a mean galactocentric luminosity distance = 338.5 Mpc, 
with the angular diameter distance thus being given by Da = -Dl/(1 + z)'^ = 288.9 Mpc. 
The luminosity and HI gas masses quoted in Table 5 were evaluated using this Dl, and the 
radial distances and scale lengths were determined using D\, with the 70" last data point 
then being at a mammoth 98.0 kpc. The rotation curve data for Malin 1 were originally 
observed in and have recently been reanalyzed in 65 1. In part because of beam smearing 
correction considerations, the authors of 65|] have revised the inner region rotation curve of 
661] quite substantially, but are in reasonable agreement with 66| in the outer region, the 
region that is of most interest to us here. For the fit to the galaxy we thus use the first four 



rotation curve data points given in 
to the 38° inclination determined by 



65| ■ and the fifth point given in [66| (as then adjusted 



65|). 



In regard to some specific galaxies within our 111 galaxy sample, we should note that 
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we have some difficulty fitting NGC 7793, with the shape of its rotation curve not readily 
lending itself to fitting. Since the HI data for this galaxy only go out to six optical disk 
scale lengths or so, the fits are very sensitive to any inner region structure that would not 
be modeled by a single exponential disk. 

For the galaxy NGC 3109 we should note that we followed [67| and scaled up the HI gas 
mass by a factor of 1.67 to allow for loss of flux in the original radio observations of the 



galaxy given in [58| . Even with this rescaling, at the one standard deviation NED distance of 
1.5 Mpc our fit still falls a little below the observed velocities at the largest radial distances. 
However, as noted in our earlier fit to this galaxy js], the fit falls right on the data at the 
slightly larger adopted distance of 1.7 Mpc. 

For the galaxy NGC 4736 the surface brightness profile was found [10] to decompose into 
a two-disk structure, a small disk with scale length 0.3 kpc that is operative in the inner 
80" = 1.9 kpc region where the ffist 13 of the 82 rotation curve data points reported in |l^ 
are located, together with a large disk with scale length 2.1 kpc scale length that is operative 
in the region greater than 80". For simplicity we opted not to truncate either of the two 
disks so that we could use equation as is for each of them. In the fitting the inner region 
disk was found to have a fitted mass 0.708 X 10^° Mq, while the dominant primary disk 
was found to have a mass 1.630 x 10^° Mq. In Table 1 the reported value for (M/L)stars 
for this galaxy is the total stellar mass of the two disks combined divided by the total blue 
luminosity of the galaxy. 

The galaxy NGC 2976 is also reported to have a two-disk structure [l^, with an effective 
-Ro = 79" = 1.4 kpc in the less than 100" radial region and an effective Rq = 34" = 0.6 
kpc in the greater than 100" radial region. With the rotation curve data of 10| ending at 
147", and with 28 of the reported 42 rotation curve data points lying in the less than 100" 
region, for the fitting we have approximated the two disks by a single disk with a blended 
scale length Rq = 1.2 kpc. 

The galaxy NGC 4826 is a highly unusual galaxy in which the inner 10 of the 89 rotation 
curve data points reported in [10[ are counter- rotating with respect to the outer 79. With 
the two regions being well segregated (the inner points lie within 50" of the center of the 
galaxy while the outer region points lie beyond 130"), we provide a fit to the 79 outer region 
points alone. 

For two of the galaxies in our sample (UGC 5999 and F571-8) we note that even though 
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the fits themselves are reasonable, we find fitted mass to light ratios much larger than the 
upper bound of 10Mq/Lq suggested by population synthesis models [55|. The galaxy UGC 
5999 has a reported inclination of 14°. However, as noted in 68|, since photometric data 
are not too sensitive to the inclination angle for close to face-on LSB galaxies, the data can 
permit a modest increase in the inclination to 22° or so. With this latter value the rotation 
velocities are reduced by a factor of sin(14°)/ sin(22°) = 0.65, with the fitted value for the 
mass to light ratio coming down to 5AM(7)/ Lq, and with the fit (not shown) even being 
improved. 

The galaxy F571-8 is close to edge on with an inclination close to 90°, and while there 
is then little sensitivity to inclination, because the galaxy is edge on, there instead are 



uncertainties in the photometry due to optical depth and projection effects [411 • For F571-8 
the rotation curve data only go out to 2.7 disk scale lengths and rise to a quite high value 
of 143.9 km sec~^. With the data being entirely in the inner optical disk region, given the 
large reported values for the velocities, a fitted disk mass as high as 4.48 x 10^*^ Mq is to be 
expected. However for this galaxy the reported luminosity is only 0.19 x 10^° Lq, to thus 
lead to a large mass to light ratio. There is thus a mismatch between the large reported 
inner region rotation velocities and the small reported luminosity. Given the photometry 
uncertainties it is possible that for this galaxy the luminosity may have been underestimated. 

Of the galaxies in the 21 large galaxy sample (the region were we are maximally sensitive 
to distance determinations) there were only three galaxies whose fitting we found challenging, 
viz. NGC 5533, NGC 6674 and UGC 2885, each a galaxy with a bulge. However, the fitting 
difficulties were mainly in the inner region where one has to make a bulge/disk decomposition 
of the luminosity, and not in the asymptotic region where the quadratic term contribution 
was still readily able to universally cancel the linear potential term contribution. Since the 
NED determination of the adopted distance for NGC 5533 is given as 47.7 ±5.7 to one 
standard deviation, we found that using 42.0 Mpc as the adopted distance gave the tightest 
fit. For NGC 6674 only a redshift-based adopted distance is available, and it lies in the 
range 42.0 to 57.0 Mpc. For this galaxy the fitting again preferred the smallest adopted 
distance value. For UGC 2885 we found that the fitting could be improved if, as described 
in the appendix, we used a bulge scale length somewhat larger than the one reported in the 



literature. For NGC 5533 and NGC 6674 we recall 



55| that NGC 5533 has significant side- 



to-side asymmetries and kinematic evidence for a warp, while NGC 6674 has a large scale 
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non-axisymmetric structure and a substantial inner region bar 45|]. Consequently we should 
not anticipate being able to do more than fit the general trend for these two galaxies in the 
inner region. Nonetheless, none of these inner region luminosity structure issues affect the 
outer region where all the various luminous components consolidate to produce one effective 
A^* in the asymptotic equation fl2^ that then readily controls the outer region. 

In Figs. 1 - 5 we present the rotational velocities with their quoted errors (in km sec~^) 
for all of the galaxies in the 111 galaxy sample as plotted as functions of radial distances 
from galactic centers (in kpc). For each galaxy we have exhibited the contribution due to 
the luminous Newtonian term alone (dashed curve), the contribution from the two linear 
terms alone (dot dashed curve), the contribution from the two linear terms and the quadratic 
terms combined (dotted curve), with the full curve showing the total contribution. As we 
see, the tightly constrained equation fl23|) captures the essence of the data, and does so 
without needing any dark matter whatsoever. 

One of the most interesting aspects of the fits is that in the galaxies that go out to 
the largest radial distances the contribution of the linear potential (dot dashed curve) would 
actually lead to an overshoot of the data, but as the Figures show this overshoot is completely 
arrested by the quadratic potential term (dotted curve). Since the quadratic term would 
eventually cause rotation velocities to fall, to illustrate the effect, for the very small DDO 154 
and for the very large UGC 128 and Malin 1 galaxies in Fig. 6 we plot the expectation of our 
model over an extended distance range. The anticipated ultimate fall in rotation velocities 
is thus a significant falsifiable diagnostic of the theory presented here, and intriguingly for 
the galaxy Malin 1 the fall is expected to set in shortly beyond the current last data point. 
For DDO 154 we note that there actually have been some suggestions of a possible fall in 
the literature. However, the small fall at the end of the rotation curve that had originally 
been reported in |69l| is not apparent in the more recent THINGS survey of the galaxy. 
Additionally, in [70| it was suggested that there might be a fall in the rotation curve at 
distances beyond those currently available. (The authors are indebted to Dr. M. Milgrom 
for alerting them to this reference.) However, the fall discussed in is thought to set in 
well before the one predicted here. 

Of particular interest in the sample are the HSB galaxies NGC 3992, NGC 3198, NGC 
2841 and UGC 2885, all four of which were also in the 21 large galaxy sample. While NGC 
3992 is part of the Ursa Major cluster study, its NED distance of 25.6 Mpc puts it well 
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beyond the 15.5 - 18.6 Mpc distance range that the Ursa Major cluster is thought to he 
within, and yet even at this much larger distance our theory is still able to accommodate 
it. Both the NGC 3198 and NGC 2841 galaxies were in the 11 galaxy sample considered in 
{sj, and the rotation curves shown here are of precisely the same shape as they had been 
previously. However, in the interim the adopted distances to both of these galaxies have 
been revised upwards by as much as 50 per cent. With the linear term contribution to v"^ 
being of the form 7oc^-R/2, it is extremely sensitive to distance determinations since 70 is 
given in equation ( 122|) as an absolute quantity. Consequently, as the Figures show, the linear 
potential terms would now be requiring the NGC 3198 and NGC 2841 rotation curves to 
rise. That no rise is seen is due entirely to the quadratic term, with the currently observed 
flatness of these rotation curves being due to a natural interplay of all the various terms 
involved. 

While the rotation curves of all of the galaxies in the sample are obtained from HI radio 
studies that extend beyond the optical disk region, the rotation curve of UGC 2885 had 
originally been obtained from HII optical studies TjJ that were thus restricted to the optical 
disk region where hot stars can ionize hydrogen gas. Now even though the UGC 2885 HII 
rotation curve data were found to quickly rise to fiat (to thereby immediately suggest a 
missing mass problem), because the optical disk is highly extended, within the optical disk 
it is actually possible to fit the UGC 2885 HII rotation curve data using only the Newtonian 
contributions of the luminous disk and bulge and visible HI gas, without the need to invoke 

n 

dark matter or alternate gravity at all [53j . Such a fit would have to be a maximum disk fit 
in which the luminous disk A^* is taken to be as large as it possibly can be in equation (I9l), 
with optical disk region flatness thus not necessarily being an indicator of any failure of the 
luminous Newtonian expectation. However, because the UGC 2885 optical disk region does 
go out so far (not relatively in disk scale lengths but absolutely in kpc), the inner region 
rotation curve is sensitive to the linear and quadratic terms in equation fl23|) . and as is seen 
in our fit to UGC 2885, they force the normalization of the Newtonian disk term to be less 
than maximal. Our work here thus supports the notion that the UGC 2885 optical disk 
region HII data do in fact serve as an indicator of the failure of the luminous Newtonian 
expectation. 

In total, our fits here and in [9] are noteworthy in that the universal 70 and k terms 
have no dependence on individual galactic properties whatsoever and yet have to work in 
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every single case. Our fits are also noteworthy in that we have captured the essence of the 
rotation curve data even though we have imposed some rather strong constraints on the input 
parameters. For adopted distances in most cases we have used NED mean values. We have 
not used actual surface brightness distributions or actual gas profiles but have treated these 
distributions simply as exponentials. Moreover, for the optical disk scale lengths we have 
mainly used those associated with the longest wavelength bands available, and have taken 
gas scale lengths to be four times disk scale lengths. Additionally, we have not included 
the effects of a disk thickness or taken any galactic inclination angle uncertainties into 
consideration. On the theoretical side our fits are noteworthy in that equation fl25]l is not 
simply a phenomenological or empirical formula that was extracted solely from consideration 
of the systematics of galactic rotation curves. Rather, equation fl25]) was explicitly derived 
from first principles in a fundamental, uniquely prescribed metric-based theory of gravity, 
namely conformal gravity. Moreover, conformal gravity itself was not even advanced for 
the purposes of addressing the dark matter problem. Rather, before it was known what its 
static, spherically symmetric solutions might even look like, it was advanced by one of us 



72| simply because it had a symmetry that could control the cosmological constant. Our 



fitting is thus quite non-trivial. 

V. GENERAL COMMENTS 

While beyond the scope of the present paper, we note that since the scale we find for k is 
of order 1/(100 Mpc)^, our work potentially has some interesting implications for clusters of 
galaxies. For clusters one can make measurements using either interior or exterior probes. 
The interior probe involves measuring galaxy kinematics and X-ray kinematics, while the 
exterior probe involves measuring lensing by clusters. For the interior case we need to use 
equation ( |T^ for points within the cluster, and for lensing we need to use equation ( fT3i) for 
points exterior to the cluster, and in both cases we need to include the global effect due to 
all of the other clusters in the Universe. Since previous applications of conformal gravity to 



clusters (velocity dispersions |73|, X-rays in clusters [7^, l75|, and lensing |76l-l80j ) did not 



include this global effect, studies of its possible impact on clusters and also on gamma ray 



bursters 



81 



484l | could be instructive. 



A second area of interest is the growth of inhomogeneities in conformal cosmology, to see 
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if one can generate a theoretical expectation for the matter distributions /(r) and h{r) that 
appear in equations (Q and f|T3|) . While a theory for inhomogeneity growth in conformal 
gravity is only in the initial stag es of development with only tensor gravitational fluctuations 
having so far been studied [85], by providing a measurement of nc^ = (1/3) /^°^ dr' r'h{r') 
in this paper we have determined one of the moment integrals of the matter source h{r). 
Our paper will thus provide an immediate test for the theory of matter fluctuations once 
it is developed. An essential first step toward developing such a theory has recently been 



taken in 



Sis, 



where it was noted that conformal cosmological perturbation theory 
has to be developed as a power series in Planck's constant rather than as a power series 
in the gravitational coupling constant. One is able to make such a quantum- mechanical 

gravity 
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development since a realization of fourth-order derivative theories such as conforma 
has recently been found in which the quantum theory is unitary and ghost free 

Another area where our theory could potentially be tested is in the behavior of satellite 
galaxies around primary galaxies. Specifically, in the above we had noted that the interplay 
of the linear and quadratic potentials in equation (l24l) would lead to a cut-off in bound 
circular orbits at the point at which R'^{d(f)/dtf = N*/3*c^/R + {N*'j + 'jo)c^R/2 - kc^R^ 
would vanish, with galaxies not being able to support bound circular orbits beyond this 
cut-off. Nonetheless, as we show in the Appendix, it is still possible for galaxies to support 
some trajectories beyond the cut-off, they just would not be circular. In a non-circular 
trajectory with energy U and angular momentum J the quantity {dr/dtY/2 + J^/2r^ = 
U + N*l3*c^/R - (iV*7 + -fo)c^R/2 + kc^R^/2 would need to be positive. Such trajectories 
are thus of relevance when the kc^R^ term is large, while the circular orbits occur when nc^R? 
is small. A possible way to explore any such large R switch over would be through the use of 
satellite galaxies as they are located well outside the primary galaxies. In making any such 
application one would have to allow for the fact that in the conformal theory the satellites 
themselves also put out linearly growing potentials, which would cause the satellite galaxies 
to interact with each other far more than they do in Newtonian gravity where potentials 
fall with distance. The requirement of a switch over in the conformal theory conforms with 
the existence of a gap between the luminous material within each primary galaxy and the 
luminous material in its satellites. Actually determining the specific way in which the switch 
over occurs could be quite instructive, and especially since there is no such switch over in 
the dark matter case, where, as noted for instance in 89|], the extrapolation for Newtonian 
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dark matter halos is quite straightforward and smooth. (We are indebted to a referee to our 
paper for alerting us to the issue of satelhte galaxies and to this reference.) 

It is also of interest to compare our work with some other alternative theories that have 
been proposed. Of them, two other non-dark-matter theories have also had success when 



applied to large samp 



theory of Milgrom 



of Moffat 
setting 



90 



es of galaxies. One is the Modified Newtonian Dynamics (MOND) 



92| , and the other is the Metric Skew Tensor Gravity (MSTG) theory 



9J]. In MOND one modifies the connection between acceleration and force by 



a \ V 



ao/ r 



dV 
dr 



(25) 



where a = v'^/r is the ordinary centripetal acceleration, dV/dr is the standard Newtonian 
gravitational force, and ^{a/ao) is the modification as defined in terms of some new universal 
parameter ao with the dimensions of acceleration. Milgrom introduced this modification 
because of his empirical discovery that in all those cases where the standard Newtonian 
theory needed dark matter, the measured centripetal accelerations were found to fall below 
a common value ag = 1.2 x 10^^ cm sec~^. Through use of the simple expression 

X 



(26) 



(1 + X2)V2 

for the function yu(a;), Milgrom was able to construct a function that interpolated between 
standard Newton-Kepler behavior at a ^ ao and departures from it in the a ^ ao MOND 
regime where it led to asymptotically flat rotation velocities. In the years since Milgrom 
first introduced MOND many rotation curves of many different varieties of galaxy have been 
measured, and to a remarkable degree (see e.g. 20|, l55|, |67|, |68|, |95|] and references therein) 
they have been successfully fitted by equation ( 126|) without the need to include any dark 
matter at all. 

In Moffat's MSTG theory a skew-symmetric tensor field is coupled to Einstein gravity, 
with the centripetal accelerations that result being given by the simple formula 

•2 GM /Mn\^/2 
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r \ r 
1 + — exp 



(27) 



for a galaxy of mass M. In applications of equation ( 1271) the parameters Mq and ro are 
found to be given by Mq = 9.60 x IO^^Mq and ro = 4.30 x 10^^ cm; and together they 
combine with Newton's constant G to give a universal acceleration parameter GMq/tq = 
6.90 X 10^^ cm sec~^. In equation (!27|l the velocity obeys v"^ = GM/r for r ^ ro and obeys 
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= GM[1 + (Mo/M)^/^]/r for r ^ tq, to thus be Kepler in both hmits, albeit with different 
effective Newton constants. Via equatio n fl2711 successful fitting to a wide variety of galaxies 
has been obtained without dark matter [96|. 

That conformal gravity, MOND and MSTG can all succeed in fitting the data is because 
not only does each one of them possess a universal (i.e. galaxy independent) parameter 
with the dimensions of an inverse length (viz. oq/c^ = 1.33 x IQ-^^ cm for MOND, 
GoMo/ry = 7.67 x 10-^9 cm^^ for MSTG, and 70 = 3.06 x IQ-^o cm^^ for conformal 
gravity), the data do too. Specifically, in the Tables we have listed the value of the quantity 
(f^/c^i?)iast at the last data point for each of the galaxies in the sample. As we see, despite 
the huge variation in luminosity and central surface brightness across the sample, within 
one order of magnitude all the (f ^/c^-R)iast values cluster around a value of 3 x 10^^*^ cm~^ 
or so. (In all the galaxies where (f ^/c^-R)iast is greater than 10 x 10~^° cm~^, the luminous 
Newtonian contribution is dominating ffast) with those galaxies not being asymptotic enough 
to be in the region where the universal linear potential 70 term would dominate.) Now 
different theories cannot agree for ever, and since equations (123|) . (!25|) and ( 127|) predict 
differing behaviors at large R, study of rotation curves at large enough R could enable us 
to distinguish between them. 

As regards the near universality of (f ^/c^i?)iast, we should note that this is an empirical 
property of the raw data themselves. Moreover, while there may be some uncertainties in 
the adopted distances to the galaxies, such uncertainties are never more than a factor of 
two or so. With the velocities being uncertain to no more than 10 to 20 per cent or so, the 
near universality of {v^/c^R)iast is thus a genuine property of the data. It should thus be 
regarded as an important empirical clue for galactic dynamics. 

It is important to recognize that the fits provided by conformal gravity (and likewise by 
MOND and MSTG) are predictions. Specifically, for all these theories the only input one 
needs is the optical and gas spectra, and the only free parameter is the M/L ratio for each 
given galaxy, with rotation velocities then being determined. Moreover, the M/ L ratios are 
highly constrained by the data in the inner rotation curve region where departures from the 
Newtonian expectation are at their minimum, and as the Tables show, they are all by and 
large found to be of order the mass to light ratio found in the local solar neighborhood, 
just as one would want. It is important to stress this point since dark matter fitting to 
galactic data works very differently. There one first needs to know the velocities so that one 
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can then ascertain the needed amount of dark matter, i.e. in its current formulation dark 
matter is only a parametrization of the velocity discrepancies that are observed and is not 
a prediction of them. Dark matter theory has yet to develop to the point where it is able to 
predict rotation velocities given a knowledge of the luminous distribution alone (or explain 
the near universality found for (f ^/c^-R)iast)- Thus dark matter theories, and in particular 
those theories that produce dark matter halos in the early Universe, are currently unable 
to make an a priori determination as to which halo is to go with which particular luminous 
matter distribution, and need to fine-tune h alo p arameters to luminous parameters galaxy by 



galaxy. (In the NFW CDM simulations 97|, |98j for instance, one finds generic spherical halo 



profiles close in form to cr(r) = o"o/['^(^+'^o)^] (as then cut off at some r^nux), but with the halo 
parameters needing to be fixed galaxy by galaxy.) No such fine-tuning shortcomings appear 
in conformal gravity, and if standard gravity is to be the correct description of gravity, then 
a universal formula akin to the one given in equation ( 123|) would need to be derived by dark 
matter theory. However, since our study establishes that global physics has an influence on 
local galactic motions, the invoking of dark matter in galaxies could potentially be nothing 
more than an attempt to describe global physics effects in purely local galactic terms. 

The authors wish to thank Dr. J. R. Brownstein, Dr. W. J. G. de Blok, Dr. J. W. Mof- 
fat, and Dr. S. S. McGaugh for helpful communications, and especially for providing their 
galactic data bases. This research has made use of the NASA/IPAC Extragalactic Database 
(NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technol- 
ogy, under contract with the National Aeronautics and Space Administration. 

Appendix A: Galaxies with Bulges or without Photometry 
1. Spherical Bulge Formalism 

For a spherically symmetric matter distribution with radial matter number density a{r) 
and N = An J dr' r''^a{r') stars, as follows directly from equation f|T2l) and equation f|T^ the 
rotational velocities associated with the Newtonian and linear potentials are given by 
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Ordinarily it is not the 3- dimensional a{r) which is directly measured in spherical astronom- 
ical systems. Rather, it is only the two-dimensional surface matter distribution I{R) which 
is measured, with a{r) having to be extracted from it via an Abel transform 



air 



IT 



dR 



r{R) 



(i?2 _r2)l/2^ 



HR) =2rdr 

I[K) ZJ^ _ ^2)1/2- 



In terms of I{R) the Newtonian integral in equation (lAip can be rewritten as 53| 
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For the very convenient exponential surface density 
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considered in 



99| . the Abel transform can be performed analytically, to yield 
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For particles orbiting such a spherical bulge at radius r we immediately obtain circular 
velocities of the form 
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2. Applications to Galaxies with Spherical Bulges 



We considered bulges for 11 galaxies. At our adopted distances the measured bulge scale 



lengths for NGC 801 and NGC 2998 are t = 0.9 kpc and t = 0.9 kpc 



55 



99|, for NGC 5371, 
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NGC 5533 and NGC 6674 the measured values are t = 0.9 kpc, t = 1.3 kpc, and t = 0.9 

n n 

kpc [55], and for ESQ 0140040 t = 1.36 kpc [42]. For NGC 5033, NGC 5055, NGC 5907 
and Malin 1 we determined respective best values of t = 1.73 kpc, t = 0.35 kpc, t = 1.84 
kpc, and t = 1.0 kpc from fitting the rotation curve themselves. The fitting could generally 
accommodate fairly broad ranges around these particular fitted values, and could do so while 
only affecting the fitting in the inner rotation curve region. (The expression for f|(r) given 
in equation f lA7p peaks at around r = 2.7t and becomes Keplerian by about r = 5t. Thus 
for all but the innermost of the points on the rotation curve, acts just like a point 

Newtonian source at the center of the galaxy. In addition, just like a linear potential point 
source, in the innermost region the contribution of the f^(r) term given in equation (lASp is 
negligible.) For UGC 2885, we could readily fit the outer 16 of the 19 rotation curve data 



points using the t = 0.6 kpc scale length given in |55|, |99|. Given the uncertainties inherent 
in bulge/disk decompositions, we can vary the bulge scale length somewhat, to find that 
we can improve the fit for the innermost three points while still being able to account for 
the other 16 points. In the Figures we report the fit with t = 1.0 kpc. In the fitting we 
obtained fitted bulge masses for the 11 galaxies that are respectively given by 4.29, 1.93, 
2.38, 11.12, 10.44, 3.52, 9.75, 0.73, 7.71, 9.46, and 8.72 (in units of IO^^Mq). While there 
may be some uncertainties in bulge/disk decompositions, these uncertainties only affect the 
inner rotation curve region and do not impact on the behavior of rotation curves at the 
largest radial distances where the missing mass problem is the most pronounced and where 
the linear and quadratic potential terms are dominant. 



3. Treatment of Galaxies with no Photometry 

For two of the galaxies in our sample (NGC 7137 and UGC 477) there appears to be 
no surface photometry reported in the literature, while for five of them (ESO 840411, ESO 
1200211, ESO 3020120, ESO 3050090 and ESO 4880490) there is only a minimal amount. In 
the absence of any surface brightness photometry our strategy is to assume that the surface 
brightness can be described as a disk with exponential S(i?) = Soe^^/^°, and simply do a 
fit to the rotation curve data using Rq and A^* as two free parameters. To constrain such fits 



we follow 



55| and require that M/L not be less than 0.2Mq/ Lq. Additionally, we require 



that Rq be less than the measured value of -Riast at the last data point. On imposing these 
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constraints, we will regard a fit as acceptable, though of course only indicative, if we can 
find a range of such constrained values of Rq and A^* for which the fitting is reasonable. 
Interestingly, this prescription is found to work for all seven of the galaxies, with there 
being a range of allowed values in each case. In the Tables we present some typical fitted 
values within the allowed ranges for each of the seven galaxies, and then use these values 
to generate the associated Figures. Just as with the spherical bulges, we should note that 
none of these photometry concerns affect the outer region rotation curve fitting. 

To support the Rq values that we obtained this way, we note that for the 
five ESO galaxies listed above, some limited surface brightness data actually are 
available. Specifically, in the ESO Lauberts-Valentijn Archive (as accessed at 



http : / / archive . eso . org / wdb / wdb / eso / esolv / form ) both a red band total apparent magni- 



tude rriT and a red band mean central surface brightness mo (in magnitudes per square 
arc second) are listed for the ESO 840411, ESO 1200211, ESO 3020120 and ESO 4880490 
galaxies, while a red band mean central surface brightness is listed for ESO 3050090. The 
quantity rho is not precisely the apparent central surface brightness itlq itself, but rather the 
average apparent surface brightness in a 10 arc second circular aperture. If we nonetheless 
now approximate rriQ by rfiQ, then from 2.51ogio(27r-Ro) ~ "^o ~ (a quantity that conve- 
niently is not affected by extinction corrections) we can extract an approximate value for Rq 
in arc seconds. Doing this is found to yield red band scale lengths i?o(ESO 840411) = 9.4", 
i?o(ESO 1200211) = 18.9", i?o(ESO 3020120) = 13.5", and i?o(ESO 4880490) = 11.1"; and 
thus respective scale lengths of 3.8, 1.4, 4.6, and 1.6 kpc at the adopted distances listed in 
the Tables. For ESO 3050090 only a blue band total apparent magnitude of 13.08 is listed. 
Taking the red band total apparent magnitude to be equal to 13.5, 13.0 and 12.5 (i.e. to be 
within 0.5 magnitudes of the blue band value, a reasonable enough expectation) respectively 
yield i?o(ESO 3050090) = 15.9", i?o(ESO 3050090) = 20.0" and /2o(ESO 3050090) = 25", 
with the 20" value corresponding to -Ro(ESO 3050090) = 1.3 kpc at the adopted distance 
listed in the Tables. For all five of the ESO galaxies then, the Ro values are compatible with 
the allowed ranges of values for Ro that we found from fitting the rotation curves. 
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4. Double-Counting in the Bulge-Disk Overlap Region 



Since bulges and disks of spiral galaxies overlap in the galactic center region, there could 
be some double counting. A possible way to allow for this would be to truncate the bulge 
contribution so that it is only non-zero in the galactic center region, and another possibility 
would be to truncate the disk contribution so that it is only non-zero outside the galactic 
center region. We describe the formalism for doing this in the case where there are both 
Newtonian and linear potentials. However, for the bulge galaxies of interest to us in this 
paper, we found that neither of the two truncation procedures had that much of an impact 
on the fits (mainly because only the innermost rotation curve points could be affected by 
the bulge/disk decomposition in the first place), and only present the formalism here for 
reference purposes. 

For bulges the most straightforward truncation is of the form 

N 



■,Qe{to-r), (A9) 



with the volume density being truncated at r = to- (We truncate a{r) rather than I{R) with 
a step function since as noted in [t^, truncating I{R) with a step function would generate 
singularities in the Abel transform.) Given the truncated equation (lAOp the orbital velocities 
are then given by 

vl(r<to) = — / dzz^Ko(z), vl(r > to) = — / dzz^Ko(z) AlO 

^ nr Jo ' nr Jo 



v^irKto) = / dzz^Ko{z) 



TT JO 



dzz'^Koiz) + 

Jo 



'*'ll)-"^M7 



Snr Jo "\ / ■ 2^^2 

v'^ir > to) = / dzz^Ko{z) - / dzz'Ko{z). (All) 

For disks the most straightforward truncation is to truncate the surface density according 
to T,{R)6{R — So) so that it starts at some minimum value Sq near to the galactic center. 
Following the procedure described in |l| for the arbitrary S(-R), we find that for a truncated 
S(i?) the Newtonian potential contribution is given by 

roc poc 

Vp{R) = -271(3*0^ / dk / dR! R!T.{R!)Jo{kR!)Jo{kR) 

Jo J so 
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27if3*c^ / dk / dR' R'E{R')Jo{kR')Jo{kR) 



+271 (3* / dk / dR' R'^{R')Jo{kR')Jo{kR) 



(A12) 



/o JO 

For an exponential disk with = Soe~^/^° (and thus a truncated number of stars 

A^TR = 27rSo(-Ro + SoRo)e~^°^^° in the R > Sq region), the first of the last two integrals in 
equation ( ]A12I) can be done analytically (and leads to the Newtonian term given in equation 
([n]))- However, the dk integration range in the second of the last two integrals in equation 
(]A12p has to be broken into two separate R < R', R > R' regions. But when sq is less than 
the positions R of the points of interest for the rotation curves (as would typically be the 
case) we only need the R> R' region, with equation (]A12p then simplifying to 



Vp{R > So) 



+ 



-7rEoP*c^R 



V2i?o 
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R\ f R 
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(A13) 



where K(x^) = Jq^"^ dy (1 — x^sin^y) is the complete elliptic integral of the first kind. 

Finally, on differentiating equation (IA13P with respect to R and recalling that 

dKjx') ^ E{x') - (1 - x')K{x') 
dx x{l — x^) 

where E(a;^) = ^q^"^ dy{l — x^sin^y)^/^ is the complete elliptic integral of the second kind, 
we then obtain for the Newtonian contribution to the orbital velocities 



vl{R > So] 



Rq 



'r 



rd^R!e-^'/'^ 
Jo 
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(A15) 
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Similarly, again following [1], for the linear potential contribution we obtain the general 
V,{R) 

JO Jso 
roo roc 

dR'R'J:{R') [{R^ + R'^)Jo{kR')Jo{kR) - 2RR' Ji{kR')Ji{kR) 

dR'R'E{R') [{R^ + R'^)Jo{kR')Jo{kR) - 2RR' Ji{kR')Ji{kR) . 

(A16) 



Thus for a truncated exponential disk we obtain 



V^{R > So) = 7rSo7*c^i?i?o 
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(A17) 



for points with R > sq. Finally, on differentiating equation ( ]A17p with respect to R and 
recalling that 



we then obtain for the linear potential contribution to the orbital velocities 
vHR > So) = nJ:ol*c'R'Roh {^]kJ ^ 



\2Rq/ \2i?o 



E(x' 



x=R'/R ' 



(A19) 



For most cases of interest the disk scale length Rq will typically be much larger than the 
truncation point sq, with the e"^'/"^" term thus being very close to one in the integration 
ranges needed for both equation (lAlSP and equation ( ]A19p . Then, when we do approximate 
the e~^'^^ term to one, the integrals in equation flAlSP and equation (IA19P can be done 
analytically, and yield the simple and convenient expressions 



vl{R > So) 
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(A20) 
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5. The Effect of Mass Generation on Galactic Orbits 

In a strictly conformal invariant theory such as conformal gravity there can be no mass 
scales at the level of the Lagrangian, with all such scales needing to be generated dynamically. 
Similarly, in the SU{3) x SU{2) x U{1) standard model of elementary particle physics there 
are no fundamental mass scales in the gauge boson and fermion sectors of the theory, with 
mass scales being generated by scalar fields that could either be fundamental q-number 
fields or c-number matrix elements of fermion bilinear condensates. Since the double-well 
scalar field potential with its tachyonic mass would violate the conformal symmetry, in the 
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conformal case the symmetry breaking must be done by fermion condensates, but in both 
the conformal and the SU{3) x SU (2) x U (1) cases there must be some form of scale-breaking 
scalars. (As noted in it is because of the contribution of the energy and momentum of 
the scale breaking fields that the traceless matter field energy-momentum tensor required 
of a conformal theory is able to support a non-zero dynamically induced matter field mass 
term.) Moreover, when one couples to grand- unified theories or to cosmology, there must be 
additional scalar fields that set the scale for the early Universe. These latter fields are not 
involved in particle mass generation (they are already operative at temperatures above the 
early Universe SU{2) x f/(l) phase transition) and can be modeled by expectation values 
of fermion quadrilinears of the form y summed over all the fermions in the 

Universe. Masses for individual fermions would be generated by bilinear condensates of 
the form (tpiipi), while masses of systems of fermions would be generated by 
as summed over the N fermions (i.e. the mass of an A^-particle system is of order 
times that of each individual particle in it). On stepping up from microscopic fermions 
to macroscopic systems such as stars and galaxies, we see that if we associate an effective 
scalar field with each solar mass unit of material, then the scalar field associated with an 
entire galaxy of mass N*Mq would be A^* times larger. Since it is this latter scalar field 
that sets the scale for the galactic geometry (by contributing to the source integrals given 
in equation ([7]) as it generates galactic masses in the first place and enables particles to 
localize), it will completely overwhelm the contribution that a one solar mass generating 
scalar field could make to the galactic orbit of a one solar mass object. (As noted in lOOj, 
the contribution of a radial scalar field S{r) to the source function /(r) in equation is 
given by /(r) = {SS" — 2S"^)/4ag, with it thus being the departure of the scalar field from 
a constant value that generates the potential produced by a localized source, with the /(r) 
of an A^-particle system being A^ times larger than the /(r) of a one-particle system.) 

To explore the expli cit effect that a scalar field has on a particle orbit we follow the 



discussion given in |lOO| and introduce the test particle action 

lT = -hJ drSix) (A22) 

that describes the coupling of a particle to its own scalar field as it propagates in a back- 
ground geometry. The utility of this action is that it not only reduces to the conventionally 
used massive test particle action when S{x) is constant, for varying S{x) it is actually fully 
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conformal invariant, since the e~"^^^ change in the scalar field is compensated by an accom- 
panying e"*^^-* change in the proper time. Variation of the action of equation (IA22I) with 
respect to the coordinates of the test particle thus yields trajectories that are conformal 
invariant, with their specific form being given as [l| 



, ^ / d'^x^ . dx^" dx'"'' 
hS I + 



dr' 



dr dr 



dr dr 



(A23) 



where the proper time dr = {—gi^^ydx^dx^Y^'^ is positive in test particle orbits. 

For the general static, spherically symmetric metric of the form rfr^ = B{r)(?dt^ — 
A{r)dr'^ — r'^dQ2, and for a scalar field that only depends on the radius r, the four equations 
of motion contained in equation (IA23P take the form 

cS' dt dr 
S dr dr^ 
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(A24) 



with the prime denoting differentiation with respect to r. Equatorial plane solutions can 
be found in which 6 is fixed to 6 = 7r/2, with the equations of motion for the three other 
coordinates given in equation flA24p then being found to admit of exact first integrals. Since 
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dr 
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d'^r ^ A' /drV ^ 
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rfr2 2A\dTj 

the structure of these first integrals depends on whether or not dr/dr is zero. In trajectories 

with non-zero dr/dr in which dcp/dr is non-zero too the integrals are of the form 

,dt 
dr 



S' + AS' 
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(A26) 
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where C, D and K are integration constants. In circular orbits in which the radius r is fixed 
the integrals are of the form 
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(A29) 
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where E and F are integration constants. In radial trajectories in which the angle is fixed 
the integrals are of the form 



cBS^ = G, 
dr 



(A30) 
(A31) 



where G and H are integration constants 

On eliminating the dependence on dr from the above sets of equations, for trajectories 
with both dr/dr and d(f)/dT non-zero we obtain 



AG'^ (dr' 



c^B^ \ dt 
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" ~b' 



AK^ dr 



B 

d(f) _ K 
cB~dI ~ 



(A32) 

(A33) 
(A34) 



with equation flA33p correcting a typographical error in lOOj]. Similarly, in circular orbits 
we obtain 

1 F 

(A35) 



ld(p _F 
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and in radial trajectories we obtain 



MP (dr\ _G^_rr 



(A36) 



In addition, we note that substituting equations flA26p . flA27D . flA28D . flA29p . flA30p and 
(1A3T|) into -gf,^{dx>'/dT){dx''/dT) = 1 yields 



D = 



in trajectories with both dr/dr and d(f)/dT non-zero, and yields 



BE^ - r^F^ 



in circular orbits. In radial trajectories we obtain 



if = 0, 



with the radial velocity then being given by 
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For the circular orbits, on using equations f lA29P and f lA38P to express E and F in terms 
of the scalar field we obtain 

1/2 



rd(f) rB' rS'{2B-rB') 



(A41) 



cdt \2 2{S + rS') J ' 

an expression for circular orbits that is exact and without approximation. Finally, with 
the one-particle 5" term being negligible in the presence of the B' that is generated by 
an A^-particle background geometry, for a one-particle system in an A^-particle background 
equation (lA4ip reduces to the standard circular orbit geodesic equation 

With the dependence on the metric coefficient A{r) dropping out of equation (1A24P in orbits 
in which dr/dr = 0, and with A{r) not appearing at all in equation (]A29p . equations (]A4ip 
and (IA42P will hold in geometries in which A{r) has any form whatsoever. In particular, 
equations flA4ip and (]A42p will hold in the geometry explicitly considered in this paper in 
which A{r) = 1/B{r). We thus justify the use of equation f lA42p in the fits presented in this 
paper. 

As regards the circular orbits, we additionally note that since equations ( 1A4ip and (]A42p 
are generic equations that hold for any radial scalar field and any choice of B{r) and A{r), 
these equations will hold for circular orbits not just in conformal gravity but in standard 
Einstein gravity as well. Thus even in standard gravity in the presence of SU (3) x SU (2) x 
f/(l) scalar fields, to order the circular orbits of one-particle systems in A^-particle 
background geometries can be well-approximated by equation ( ]A42p . The only situation in 
which one could not immediately drop the S"-dependent term in f lA4ip would be in a few- 
particle system, and it would thus be of interest to see if binary pulsar systematics might 
show some sensitivity to the scalar fields that are responsible for the masses that the two 
stars in the binary possess. 

To recover the standard non-relativistic limits for the above orbits and trajectories, to 
first order in weak gravity we set 5^/(7^ = 1 — 2f//c^, J = cK/C, and B{r) = 1 + 2V{r)/c^ 
where U, J and V{r) are the first-order energy, angular momentum and potential of the test 
particle. On requiring only that A{r) = 1 + 0{V{r)) (in both the conformal gravity solution 
and in the standard Schwarzschild solution we have A{r) = 1/ B{r) ~ 1 — 2V"(r)/c^), then 
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from equations (1A32P and (1A34P as evaluated with D = we obtain 



to lowest order for a unit mass particle in a non-circular trajectory. On differentiating 
equation (IA43P we find that in a non-circular trajectory the acceleration is given by 



As regards weak gravity circular orbits, we note that since the source function /(r) in 
equation ([7]) is given by /(r) = {SS" — 2S"^)/4ag, the departure of B{r) from its flat 
spacetime value of one is correlated with the departure of S{r) from a constant value. With 
the leading term in S{r) being the mass of the test particle, for weak gravity the rS'/S term 
in equation (IA41I) is negligible, with equation (IA41I) then reducing to 

r^f^V = r^ (A45) 
\dt J dr 

for circular orbits in weak gravity (to thus also satisfy equation flA44p since d'^r/dt'^ = 0). 

For the generic potential V{r) = — /3c^/r + 7c^r/2 — /tc^r^/2 of interest to us in this paper, 
we see that for circular orbits we require 

^ + 2^ _ ^cV2 > 0, (A46) 
r 2 

while for non-circular trajectories we require 

t/ + ^- I^ + ^>0. (A47) 

Thus when the nc^r"^ term is small the conformal theory supports circular orbits, and when 
the KC^r"^ term is large the theory supports non-circular trajectories, with the non-circular 
trajectories not specifically needing to be purely radial. 

Added note. Since completing this paper we have extended our study of rotation curve 
fitting based on equatio n (1231) to include 27 additional galaxies, the bulk of which are dwarf 



spirals. As reported in |101| . on using the same values for 7*, 70 and k that are used here, 
and again with no need for any dark matter whatsoever, we find acceptable fits of the same 
quality as the ones we present here. This latest study brings to 138 the number of rotation 
curves of galaxies that have successfully been fitted by the conformal gravity theory. 
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TABLE I: Properties of the THINGS 18 Galaxy Sample 



Galaxy 


Type 


Distance 


Lb 


Ro 


-Rlast 


Mhi 




(M/L),tars 


{vyc^R)i^t 


Data Sources 






(Mpc) 


(lOiOLo) 


(kpc) (kpc) (IO^OMq) 


(lOiOAfo) 


(Mo/Lo) 


(10^30 cjjj-l) 


V 


L 


Ro HI 


DDO 0154 


LSB 


4.2 


0.007 


0.8 


8.1 


0.03 


0.003 


0.45 


1.12 


m 




[12] M 


IC 2574 


LSB 


4.5 


0.345 


4.2 


13.1 


0.19 


0.098 


0.28 


1.69 


[10] 


m 


[13] [llj 


NGC 0925 


LSB 


8.7 


1.444 


3.9 


12.4 


0.41 


1.372 


0.95 


4.17 


[10] 




[12] M 


NGC 2403 


HSB 


4.3 


1.647 


2.7 


23.9 


0.46 


2.370 


1.44 


2.89 


[10] 


\m 


M M 


NGC 2841 


HSB 


14.1 


4.742 


3.5 


51.6 


0.86 


19.552 


4.12 


5.83 


[10] 


\m 


[15] M 


NGC 2903 


HSB 


9.4 


4.088 


3.0 


30.9 


0.49 


7.155 


1.75 


3.75 


[10] 


\m 


M M 


NGC 2976 


LSB 


3.6 


0.201 


1.2 


2.6 


0.01 


0.322 


1.60 


10.43 


[10] 


M 


[16] M 


NGC 3031 


HSB 


3.7 


3.187 


2.6 


15.0 


0.38 


8.662 


2.72 


9.31 


[10] 






NGC 3198 


HSB 


14.1 


3.241 


4.0 


38.6 


1.06 


3.644 


1.12 


2.09 


m 


\m 


M M 


NGC 3521 


HSB 


12.2 


4.769 


3.3 


35.3 


1.03 


9.245 


1.94 


4.21 


[10] 


\m 


[12] M 


NGC 3621 


HSB 


7.4 


2.048 


2.9 


28.7 


0.89 


2.891 


1.41 


3.18 


[10] 


\m 


[10] M 


NGC 3627 


HSB 


10.2 


3.700 


3.1 


8.2 


0.10 


6.622 


1.79 


15.64 


[10] 


M 


[12] M 


NGC 4736 


HSB 


5.0 


1.460 


2.1 


10.3 


0.05 


1.630 


1.60 


4.66 


[10] 




[10] M 


NGC 4826 


HSB 


5.4 


1.441 


2.6 


15.8 


0.03 


3.640 


2.53 


5.46 


[10] 


\m 


[18] M 


NGC 5055 


HSB 


9.2 


3.622 


2.9 


44.4 


0.76 


6.035 


1.87 


2.36 


[10] 




[12] \U] 


NGC 6946 


HSB 


6.9 


3.732 


2.9 


22.4 


0.57 


6.272 


1.68 


6.39 


[10] 




[12] \U] 


NGC 7331 


HSB 


14.2 


6.773 


3.2 


24.4 


0.85 


12.086 


1.78 


9.61 


[10] 


[yj 


[12] M 


NGC 7793 


HSB 


5.2 


0.910 


1.7 


10.3 


0.16 


0.793 


0.87 


3.61 


[10] 




[12] M 
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TABLE II: Properties of the Ursa Major 30 Galaxy Sample 



Galaxy 


Type 


Distance 


Lb 


Rq 


-Rlast 


Mm 


-^disk 


(M/L),tars 


{vyc^R)i^t 


Data Sources 






(Mpc) 


(lOiOLo) 


(kpc) (kpc) (IO^OMo) 


(lOiOMo) 


(Mo/Lo) 


(10-30 cm-i) 


V 


L 


Rq 


HI 


NGC 3726 


HSB 


17.4 


3.340 


3.2 


31.5 


0.60 


3.82 


1.15 


3.19 


[19] 


[20] 




[20] 


NGC 3769 


HSB 


15.5 


0.684 


1.5 


32.2 


0.41 


1.36 


1.99 


1.43 


[19J 


m 


m 


m 


NGC 3877 


HSB 


15.5 


1.948 


2.4 


9.8 


0.11 


3.44 


1.76 


10.51 


Ll9j 


[20] 




[20] 


NGC 3893 


HSB 


18.1 


2.928 


2.4 


20.5 


0.59 


5.00 


1.71 


3.85 


[19] 


[20] 




[20] 


NGC 3917 


LSB 


16.9 


1.334 


2.8 


13.9 


0.17 


2.23 


1.67 


4.85 


[19] 


[20] 




[20] 


NGC 3949 


HSB 


18.4 


2.327 


1.7 


7.2 


0.35 


2.37 


1.02 


14.23 


[19] 


[20] 




[20] 


NGC 3953 


HSB 


18.7 


4.236 


3.9 


16.3 


0.31 


9.79 


2.31 


10.20 


[19] 


[20] 




[20] 


NGC 3972 


HSB 


18.6 


0.978 


2.0 


9.0 


0.13 


1.49 


1.53 


7.18 


[19] 


m 




m 


NGC 3992 


HSB 


25.6 


8.456 


5.7 


49.6 


1.94 


13.94 


1.65 


4.08 


[19] 


[20] 




[20] 


NGC 4010 


LSB 


18.4 


0.883 


3.4 


10.6 


0.29 


2.03 


2.30 


5.03 


[19] 


[20] 




[20] 


NGC 4013 


HSB 


18.6 


2.088 


2.1 


33.1 


0.32 


5.58 


2.67 


3.14 


[19] 


[20] 




[20] 


NGC 4051 


HSB 


14.6 


2.281 


2.3 


9.9 


0.18 


3.17 


1.39 


8.52 


[19] 


[20] 


[22] 


[20] 


NGC 4085 


HSB 


19.0 


1.212 


1.6 


6.5 


0.15 


1.34 


1.11 


10.21 


[19] 


[20] 




[20] 


NGC 4088 


HSB 


15.8 


2.957 


2.8 


18.8 


0.64 


4.67 


1.58 


5.79 


[19] 


[20] 




[20] 


NGC 4100 


HSB 


21.4 


3.388 


2.9 


27.1 


0.44 


5.74 


1.69 


3.35 


[19] 


m 




m 


NGC 4138 


LSB 


15.6 


0.827 


1.2 


16.1 


0.11 


2.97 


3.59 


5.04 


[19] 


m 




m 


NGC 4157 


HSB 


18.7 


2.901 


2.6 


30.9 


0.88 


5.83 


2.01 


3.99 


[19] 


[20] 




[20] 


NGC 4183 


HSB 


16.7 


1.042 


2.9 


19.5 


0.30 


1.43 


1.38 


2.36 


[19] 


[20] 




[20] 


NGC 4217 


HSB 


19.6 


3.031 


3.1 


18.2 


0.30 


5.53 


1.83 


6.28 


[19] 


[20] 




[20] 


NGC 4389 


HSB 


15.5 


0.610 


1.2 


4.6 


0.04 


0.42 


0.68 


9.49 


[19] 


[20] 


m 


[20] 


UGC 6399 


LSB 


18.7 


0.291 


2.4 


8.1 


0.07 


0.59 


2.04 


3.42 


[19] 


[20] 




[20] 


UGC 6446 


LSB 


15.9 


0.263 


1.9 


13.6 


0.24 


0.36 


1.36 


1.70 


[19] 


[20] 


[22] 


[20] 


UGC 6667 


LSB 


19.8 


0.422 


3.1 


8.6 


0.10 


0.71 


1.67 


3.09 


[19] 


[20] 


m 


[20] 


UGC 6818 


LSB 


21.7 


0.352 


2.1 


8.4 


0.16 


0.11 


0.33 


2.35 


[19] 


[20] 


m 


[20] 


UGC 6917 


LSB 


18.9 


0.563 


2.9 


10.9 


0.22 


1.24 


2.20 


4.05 


[19] 


[20] 




[20] 


UGC 6923 


LSB 


18.0 


0.297 


1.5 


5.3 


0.08 


0.35 


1.18 


4.43 


[19] 


[20] 


[22] 


[20] 


UGC 6930 


LSB 


17.0 


0.601 


2.2 


15.7 


0.29 


1.02 


1.69 


2.68 


[19] 


[20] 




[20] 


UGC 6973 


HSB 


25.3 


1.647 


2.2 


11.0 


0.35 


3.99 


2.42 


10.58 


[19] 


[20] 


[22] 


[20] 


UGC 6983 


LSB 


20.2 


0.577 


2.9 


4ff.6 


0.37 


1.28 


2.22 


2.43 


[19] 


[20] 




[20] 


UGC 7089 


LSB 


13.9 


0.352 


2.3 


7.1 


0.07 


0.35 


0.98 


3.18 


[19] 


[20] 




[20] 



TABLE III: Properties of the LSB 20 Galaxy Sample 



Galaxy 


Type 


Distance 


Lb 


Rq 


-Rlast 


AfHI 


Mdisk 


(M/L)stars 




Data S 


ources 






(Mpc) 


(IOIOLq) 


(kpc) (kpc) 


(lOiOAfo) 


(lOiOMo) 


(Mq/Lq) 


(10-30 cm-i) 


V 


L 


Rq 


HI 


DDO 0064 


LSB 


6.8 


0.015 


1.3 


2.1 


0.02 


0.04 


2.87 


6.05 


[23] 


m 


m 


[25] 


F563-1 


LSB 


46.8 


0.140 


2.9 


18.2 


0.29 


1.35 


9.65 


2.44 


[26] 


m 


m 


m 


F563-V2 


LSB 


57.8 


0.266 


2.0 


6.3 


0.20 


0.60 


2.26 


6.15 


[29] 






[28] 


F568-3 


LSB 


80.0 


0.351 


4.2 


11.6 


0.30 


1.20 


3.43 


3.16 


[26] 






[28] 


F583-1 


LSB 


32.4 


0.064 


1.6 


14.1 


0.18 


0.15 


2.32 


1.92 


[26] 






[28] 


F583-4 


LSB 


50.8 


0.096 


2.8 


7.0 


0.06 


0.31 


3.25 


2.52 


[26] 






[28] 


NGC 0959 


LSB 


13.5 


0.333 


1.3 


2.9 


0.05 


0.37 


1.11 


7.43 


[23] 


[30] 




[30] 


NGC 4395 


LSB 


4.1 


0.374 


2.7 


0.9 


0.13 


0.83 


2.21 


2.29 


[29] 


[32] 


m 


[32] 


NGC 7137 


LSB 


25.0 


0.959 


1.7 


3.6 


0.10 


0.27 


0.28 


3.91 


[23] 


[33] 


ES 


[33] 


UGC 0128 


LSB 


64.6 


0.597 


6.9 


54.8 


0.73 


2.75 


4.60 


1.03 


m 




[35] 


[35] 


UGC 0191 


LSB 


15.9 


0.129 


1.7 


2.2 


0.26 


0.49 


3.81 


15.48 


[23] 


[36] 


[36] 


[32 


UGC 0477 


LSB 


35.8 


0.871 


3.5 


10.2 


1.02 


1.00 


1.14 


4.42 


[29] 


[30] 


ES 


[30] 


UGC 1230 


LSB 


54.1 


0.366 


4.7 


37.1 


0.65 


0.67 


1.82 


0.97 






[35] 


[35] 


UGC 1281 


LSB 


5.1 


0.017 


1.6 


1.7 


0.03 


0.01 


0.53 


3.02 


[29] 


[36] 


m 


[32] 


UGC 1551 


LSB 


35.6 


0.780 


4.2 


6.6 


0.44 


0.16 


0.20 


3.69 


[23] 


[38] 


[39] 


[38] 


UGC 4325 


LSB 


11.9 


0.373 


1.9 


3.4 


0.10 


0.40 


1.08 


7.39 


[23] 


[32] 


m 


[32] 


UGC 5005 


LSB 


51.4 


0.200 


4.6 


27.7 


0.28 


1.02 


5.11 


1.30 






[35] 


[35] 


UGC 5750 


LSB 


56.1 


0.472 


3.3 


8.6 


0.10 


0.10 


0.21 


1.58 


[29] 




[35] 


[35] 


UGC 5999 


LSB 


44.9 


0.170 


4.4 


15.0 


0.18 


3.36 


19.81 


5.79 


m 




[35] 


[35] 


UGC 11820 


LSB 


17.1 


0.169 


3.6 


3.7 


0.40 


1.68 


9.95 


8.44 


[23] 


[371 


[401 


[371 
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TABLE IV: Properties of the LSB 21 Galaxy Sample 



Galaxy 


Type 


Distance 


Lb 


Rq 


Rlast 


Mhi 




(M/L)stars 


{vyc^R)i^t 


Data S 


lources 






(Mpc) 


(lO^^Lo) 


(kpc) 


(kpc) (IO^^Mq) 


(lO^OMo) 


(Mq/Lq) 


(10-30 cm-i) 


V 


L 


Rq 


HI 


ESQ 0140040 


LSB 


217.8 


7.169 


10.1 


30.0 




20.70 


3.38 


8.29 


[26] 


m 


m 


NA 


ESQ 0840411 


LSB 


82.4 


0.287 


3.5 


9.1 




0.06 


0.21 


1.49 


[26] 


m 


ES 


NA 


ESQ 1200211 


LSB 


15.2 


0.028 


2.0 


3.5 




0.01 


0.20 


0.66 


[26] 


m 


ES 


NA 


ESQ 1870510 


LSB 


16.8 


0.054 


2.1 


2.8 




0.09 


1.62 


2.02 


[26] 


m 


[43] 


NA 


ESQ 2060140 


LSB 


59.6 


0.735 


5.1 


11.6 




3.51 


4.78 


4.34 


[26] 


m 


[42] 


NA 


ESQ 3020120 


LSB 


70.9 


0.717 


3.4 


11.2 




0.77 


1.07 


2.37 


[26] 


m 


ES 


NA 


ESQ 3050090 


LSB 


13.2 


0.186 


1.3 


5.6 




0.06 


0.32 


1.87 


[26] 


m 


ES 


NA 


ESQ 4250180 


LSB 


88.3 


2.600 


7.3 


14.6 




4.79 


1.84 


5.17 


[26] 


m 


[42] 


NA 


ESQ 4880490 


LSB 


28.7 


0.139 


1.6 


7.8 




0.43 


3.07 


4.34 


[26] 


m 


ES 


NA 


F571-8 


LSB 


50.3 


0.191 


5.4 


14.6 


0.16 


4.48 


23.49 


5.10 


[26] 


[28] 




[28] 


F579-V1 


LSB 


86.9 


0.557 


5.2 


14.7 


0.21 


3.33 


5.98 


3.18 


[26] 


[28] 




[28] 


F730-V1 


LSB 


148.3 


0.756 


5.8 


12.2 




5.95 


7.87 


6.22 


[26] 


[40] 


[40] 


NA 


UGC 04115 


LSB 


5.5 


0.004 


0.3 


1.7 




0.01 


0.97 


3.42 


[26] 


m 


m 


NA 


UGC 06614 


LSB 


86.2 


2.109 


8.2 


62.7 


2.07 


9.70 


4.60 


2.39 


[26] 


m 


[35] 


[35] 


UGC 11454 


LSB 


93.9 


0.456 


3.4 


12.3 




3.15 


6.90 


6.79 


[26] 


m 


[40] 


NA 


UGC 11557 


LSB 


23.7 


1.806 


3.0 


6.7 


0.25 


0.37 


0.20 


3.49 


[26] 


m 


[32] 


[32] 


UGC 11583 


LSB 


7.1 


0.012 


0.7 


2.1 




0.01 


0.96 


2.15 


[26] 


m 


[40] 


NA 


UGC 11616 


LSB 


74.9 


2.159 


3.1 


9.8 




2.43 


1.13 


7.49 


[26] 


m 


[40] 


NA 


UGC 11648 


LSB 


49.0 


4.073 


4.0 


13.0 




2.57 


0.63 


5.79 


[26] 


m 


[40] 


NA 


UGC 11748 


LSB 


75.3 


23.930 


2.6 


21.6 




9.67 


0.40 


1.01 


[26] 


m 


[40] 


NA 


UGC 11819 


LSB 


61.5 


2.155 


4.7 


11.9 




4.83 


2.24 


7.03 


[26] 


m 


m 


NA 



42 



TABLE V: Properties of the Miscellaneous 22 Galaxy Sample 



Galaxy 


Type 


Distance 


Lb 


Rq 


Rlast 


Mhi 


Mdisk 


(M/L)stars 




Data Sources 






(Mpc) 


(lO^OLo 


) (kpc) (kpc) 


(lOiOMo) 


(lOiOMo) 


(Mq/Lq) 


(10-30 cm-i) 


V 


L 


Ro HI 


DDO 0168 


LSB 


4.5 


0.032 


1.2 


4.4 


0.03 


0.06 


2.03 


2.22 


[45j 


[45] 


[45] [45] 


DDO 0170 


LSB 


16.6 


0.023 


1.9 


13.3 


0.09 


0.05 


1.97 


1.18 


m 


[46] 


[46] [46] 


M 0033 


HSB 


0.9 


0.850 


2.5 


8.9 


0.11 


1.13 


1.33 


4.62 


m 


m 


[48] [42 


NGC 0055 


LSB 


1.9 


0.588 


1.9 


12.2 


0.13 


0.30 


0.50 


2.22 


[49] 


[49] 


[49] [49] 


NGC 0247 


LSB 


3.6 


0.512 


4.2 


14.3 


0.16 


1.25 


2.43 


2.94 


[50] 


[51J 


[51J [50] 


NGC 0300 


LSB 


2.0 


0.271 


2.1 


11.7 


0.08 


0.65 


2.41 


2.69 


[52] 


[51J 


m [52] 


NGC 0801 


HSB 


63.0 


4.746 


9.5 


46.7 


1.39 


6.93 


2.37 


3.59 


[45] 


[45] 


[53] [45] 


NGC 1003 


LSB 


11.8 


1.480 


1.9 


31.2 


0.63 


0.66 


0.45 


1.53 


m 


m 


m [55j 


NGC 1560 


LSB 


3.7 


0.053 


1.6 


10.3 


0.12 


0.17 


3.16 


2.16 


[45] 


[45] 


[45] [45] 


NGC 2683 


HSB 


10.2 


1.882 


2.4 


36.0 


0.15 


6.03 


3.20 


2.28 


[56] 


[56] 


mi [55] 


NGC 2998 


HSB 


59.3 


5.186 


4.8 


41.1 


1.78 


7.16 


1.75 


3.43 


[45] 


[45] 


[53] [45] 


NGC 3109 


LSB 


1.5 


0.064 


1.3 


7.1 


0.06 


0.02 


0.35 


2.29 


[58] 


[59] 


[59] [58] 


NGC 5033 


HSB 


15.3 


3.058 


7.5 


45.6 


1.07 


0.27 


3.28 


3.16 


[15] 


[15] 


M [15j 


NGC 5371 


HSB 


35.3 


7.593 


4.4 


41.0 


0.89 


8.52 


1.44 


3.98 


[15] 


[15] 


M [15] 


NGC 5533 


HSB 


42.0 


3.173 


7.4 


56.0 


1.39 


2.00 


4.14 


3.31 


[45] 


[45] 


m m 


NGC 5585 


HSB 


9.0 


0.333 


2.0 


14.0 


0.28 


0.36 


1.09 


2.06 


m 


m 


[60] [60] 


NGC 5907 


HSB 


16.5 


5.400 


5.5 


48.0 


1.90 


2.49 


1.89 


3.44 


m 


[55] 


[61] [55] 


NGC 6503 


HSB 


5.5 


0.417 


1.6 


20.7 


0.14 


1.53 


3.66 


2.30 


[15] 


[15] 


M [15] 


NGC 6674 


HSB 


42.0 


4.935 


7.1 


59.1 


2.18 


2.00 


2.52 


3.57 


[45] 


[45] 


M [45] 


UGC 2259 


LSB 


10.0 


0.110 


1.4 


7.8 


0.04 


0.47 


4.23 


3.76 


[62] 


[62] 


[63] [62] 


UGC 2885 


HSB 


80.4 


23.955 


13.3 


74.1 


3.98 


8.47 


0.72 


4.31 


m 


[53] 


[53] [55] 


Malin 1 


LSB 


338.5 


7.912 


84.2 


98.0 


5.40 


1.00 


1.32 


1.77 


[65] 


[65] 


[66] [65] 
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FIG. 1: Fitting to the rotational velocities (in km sec~^) of the THINGS 18 galaxy sample with 
their quoted errors as plotted as a function of radial distance (in kpc). For each galaxy we have 

exhibited the contribution due to the luminous Newtonian term alone (dashed curve), the 
contribution from the two linear terms alone (dot dashed curve), the contribution from the two 
linear terms and the quadratic terms combined (dotted curve), with the full curve showing the 

total contribution. No dark matter is assumed. 
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FIG. 2: Fitting to the rotational velocities of the Ursa Major 30 galaxy sample - Part 1 
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FIG. 2: Fitting to the rotational velocities of the Ursa Major 30 galaxy sample - Part 2 
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FIG. 3: Fitting to the rotational velocities of the LSB 20 galaxy sample 
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FIG. 4: Fitting to the rotational velocities of the LSB 21 galaxy sample 
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FIG. 5: Fitting to the rotational velocities of the miscellaneous 22 galaxy sample 
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FIG. 6: Extended distance predictions for DDO 154, UGC 128, and Malin 1. 
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